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ABSTRACT 

The Next Generation Virgo Cluster Survey has recently determined the luminosity function of 
galaxies in the core of the Virgo cluster down to unprecedented magnitude and surface brightness 
limits. Comparing simulations of cluster formation to the derived central stellar mass function, we 
attempt to estimate the stellar-to-halo-mass ratio (SHMR) for dwarf galaxies, as it would have been 
before they fell into the cluster. This approach ignores several details and complications, e.g., the 
contribution of ongoing star formation to the present-day stellar mass of cluster members, and the 
effects of adiabatic contraction and/or violent feedback on the subhalo and cluster potentials. The final 
results are startlingly simple, however; we find that the trends in the SHMR determined previously 
for bright galaxies appear to extend down in a scale-invariant way to the faintest objects detected 
in the survey. These results extend measurements of the formation efficiency of field galaxies by two 
decades in halo mass or five decades in stellar mass, down to some of the least massive dwarf galaxies 
known, with stellar masses of ^ IO^Mq. 

Subject headings: dark matter - galaxies: clusters: individual (Virgo) ~ galaxies: dwarf -- galaxies: 

formation - galaxies: general - galaxies: luminosity function, mass function 


1. INTRODUCTION 

Galaxy formation is complicated; even a sketch de¬ 
scription of the process involves cosmological structure 
formation in the dark matter component, gas cooling, 
disk formation, molecular cloud formation, star forma¬ 
tion and stellar feedback from winds and supernovae, 
central black hole growth, feedback from an active nu¬ 
cleus, secular evolution, environmental effects, mergers, 
cold gas accretion, and the contributio n from an inter- 
galactic ionizing background (see e.g.. iMo et al.l[2?Hfll 
for an overview and further references). Each of these 
elements in turn involves complex, multi-scale physics 
that we are only starting to appreciate. Faced with this 
complexity, ab initio models of galaxy formation must 
necessarily be calibrated by comparing their global pre- 
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dictions with observations. Theoretical arguments and 
dissipationless simulations do provide a simple, robust 
framework for understanding the abundance and clus¬ 
tering of cosmological structures, however - the conven- 
tional picture of cold d ark matter (CDM) halo formation 
(|Frenk fc Whitel 120121 ). How galaxies occupy CDM ha¬ 
los is then an empirical question that can be addressed 
directly by large surveys, without explicit reference to 
the complex physics responsible for the final result. 

On intermediate mass scales, the relationship between 
galaxies and CDM halos appears to be straightforward 
and one-to-one; a galaxy like the Milky Way includes 
within itself most of the stellar mass contained in its 
surrounding dark matter halo, and any halo as mas¬ 
sive as that of the Milky Way probably hosts a similar 
dominant luminous galaxy. At the larger and smaller 
ends of the halo mass scale, however, the process of 
galaxy formation becomes more complicated and less ef¬ 
ficient. The largest (galaxy cluster) halos contain thou¬ 
sands of individual galaxies rather than a single domi¬ 
nant one, and most of their baryonic mass is in the form 
of hot gas, not stars. The smallest halos, on the other 
hand, do not appear to contain any galaxies or stars at 
all. It remains unclear exactly where in the mass spec¬ 
trum the efficiency of galaxy formation drops to zero, 
and exactly which processes suppress star formation on 
which scales. Over some intermediate range of masses, 
however, we can take the ratio of the stellar mass of 
the central galaxy to the total mass of the surrounding 
halo as a simple indicator of the net efficiency of galaxy 
formatio n. This stella r -to-ha l o-mass ratio ( S HMR 
- e.g., iLeauthaiid et al.l 120121 : iBehroozi et al.l 120131 
iMoster et al.ll2013i and references therein) has now been 
measured for a wide range of systems, either individually, 
using internal kinematics to estimate total mass (e.g., 
iBlanton et al.l[200^ iMiller et al.ll2014ll or as an average 
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in well-defined samples , using galaxy-galaxy lensing (e.g., 
Leauthaud et al.l l2012t iVelander et al.l 12014 iHan et alJ 


2015t [Hudson et alJ 20151; Shan et alj 20151). satellite 


kinematics (e.g., Conrov et al ■I 120071: IMnre et a, 1.1 
IWoitak fc Mamp nll2013D . and overal l abundance and/or 
clustering fe.g.. iBehrooz et al.ll201^ iMoster et al.ll2013l : 
iRodriguez-Puebla et alJ 1201511 to estimate the average 
halo mass for the sample. The observationally derived 
SHMR provides a solid point of comparison for mod¬ 
els of galaxy formation, bridging the gap between large 
statistical samples from surveys and detailed models or 
simulations of individual galaxies. 

From the measurements to date, the SHMR has sev¬ 
eral interesting features. It reaches a maximum, corre¬ 
sponding to a peak in the efficiency of galaxy formation, 
at a characteristic halo mass comparable to the Milky 
Way’s, Mh At higher and lower masses, 

the ratio is roughly a power law (i.e., M^:/Mh oc {Mh)^, 
where the slope a has different values for small and large 
halo masses), although at either end of the scale it re¬ 
mains poorly constrained. Observations are beginning to 
probe the evolution of the ratio with redshift; the results 
are complex, with e vidence for non-monotonic, mass- 
depen dent evolution (jBehroozi et al.l[^01,'ll : IMoster et al.l 
l2013f) . The SHMR is generally very hard to measure, par¬ 
ticularly for low-mass galaxies; at low masses and high 
redshifts the measurement is essentially impossible, since 
we cannot measure dynamics for individual systems and 
cannot get enough signal from stacked samples of dwarf 
galaxies to determine an average halo mass from lensing. 
This is regrettable, given the suggestion that the SHMR 


changes signific antly at the low-mass end as one goes t o 
higher redshift (|Beiiroozi et al.l[2013l : iMiller et al.ll2014ll . 

In the local universe, the smallest dwarf galaxies are 
highly cl ustered, occurring mainly a s satellites of brighter 
systems (jKarachentsev et al.l 120131 1. Even in the Lo¬ 
cal Group, where the dwarf population is best stud¬ 
ied, the satellites of the Milky Way and Andromeda 
seem very strongly clustered comp ared to the predicted 
distribution of ha lo substructure (jKravtsov et al.l [2004 
iTavlor et al.ll200'3 l. In the context of structure forma¬ 
tion, this biasing suggests that dwarfs are an old popu¬ 
lation, having formed independently and then been ac¬ 
creted into larger systems at high redshift. The best 
place to look for ancient, low-mass galaxies in the lo¬ 
cal universe should therefore be in the densest regions. 
The centers of the two main halos comprising the Lo¬ 
cal Group host the faintest known dwarf galaxies, but 
current samples are limited and incomplete; for larger 
and possibly older samples, we must search in the cores 
of nearby groups and clusters. The nearest large clus¬ 
ter, Virgo, contains n early 900 galaxies br ighter than 
Lb = 3 X IO^Lr (iBingEfeli et ahl 1198511 . versus 10 
in the Local Group (jMcConnachiell2ni2l) , a nd has a to¬ 
tal ii r ass of vit - p = 5.2 x IO^^Mq (jMcLaughlinI 

119991 : iFerrarese et al. 1201211 . versus M 2 ooc,lg = 2.66 x 
IO^^Mq (jvan der Marel et'mi2012tl . so it represents a 
much richer hunting ground for low-mass structure. 

The Next Generation Virgo Gluster Survey (NGVS 
- IFerrarese et al.112012L Paper I hereafter) is a wide- 
area, multi-band imaging survey of the Virgo cluster 
using MegaCam on the Ganada-France-Hawaii Tele¬ 
scope (GFHT). It represents the first major revision to 


our picture of Virgo in the optical since the work of 
Binggeli, Sandage, an d Tammann more than 25 years ago 
(|Binggeli et al.lll9^ . Other papers in the NGVS series 
related to the topics considered here includ e those on the 
distri bution of globular clusters in Virgo (jPurrell et al.1 
l2014h . the propertie s of star clusters, UGDs and galaxies 
in the cluster core (iZhu et al.1120141 : IZhang et al.1120151 : 
Liu et al. in prepa ration), the dynam ical properties of 
low-mass galaxies (|Guerou et al.l l20H . Toloba et al. in 
pre paration), interactio ns within possible infalling galax¬ 
ies (iPaudel e t al.ll2013ll and optical-IR source classifica¬ 
tion ( Munoz et al.ll2014[ l. 

Using newly developed techniq ues for flat-fieldin g and 
scattered light removal (Paper 1; iDuc et all 120151 : Guil- 
landre et al. 2015, in preparation), the survey is ex¬ 
tremely sensitive to low-mass and low-surface-brightness 
dwarfs. In principle, if we could relate the resulting stel¬ 
lar mass function to a predicted subhalo mass function, 
derived from a simulation of cluster formation, we could 
determine the SHMR at extremely low masses. Further¬ 
more, since most of the galaxies in the core of Virgo were 
incorporated into this structure long ago, they should re¬ 
flect the SHMR of held galaxies at high redshift. 

There are a number of complications in carrying out 
this comparison, however. First, galaxies in clusters cor¬ 
respond to subhalos within larger halos. To determine 
the held SHMR we need to reconstruct the original mass 
a given subhalo had when it was hrst incorporated into 
the cluster. Second, baryons in the form of gas or stars 
may affect the dynamical evolution of subhalos, helping 
them survive the tidal mass loss and disruption seen in 
GDM-only simulations. Ongoing star formation will also 
add to the stellar mass of the galaxy, so the present-day 
stellar mass may not rehect its original mass on infall. 
Finally, obtaining a reliable estimate of the CDM sub¬ 
structure mass function in the core of a dense cluster is 
numerically challenging, since these regions are subject 
to the strongest resolution effects. 

In this paper we generate a set of high-resolution simu¬ 
lations of GDM halos with properties similar to the Virgo 
cluster. From the simulations we determine the mean 
subhalo ‘infair mass function (SIMF), that is the abun¬ 
dance of subhalos as a function of the mass they had at 
infall, for subhalos now at the center of the cluster, as 
well as the distribution of ‘infall redshifts’ at which these 
objects were first incorporated into the larger structure. 
(The SIMF corre sponds to the ‘subhalo initial mass func¬ 
tion’ discussed in iTavlor fc Babulll2(104l or the ‘unevolved 
subh a lo mass function’ di s cussed in Ivan den Bosch et al. 
2005 : iGiocoli et al.112008L 120101 : iJiang fc van den Bosch 
2014f) . We introduce several different models of the SIMF 
in an attempt to correct both for resolution effects and 
for the physical differences between subhalos and galax¬ 
ies. Assuming a monotonic relationship between the pre¬ 
dicted infall mass distribution and the observed stellar 
mass distribution, we calculate the SHMR for these sys¬ 
tems when they first formed in the field. Our method 
makes a number of simplifying assumptions, some of 
which require further validation through more detailed 
simulations, but it provides an initial estimate of the ef¬ 
ficiency of galaxy formation at very low mass and mod¬ 
erate redshift. 

The outline of the paper is as follows. In section [2] we 
summarize the NGVS and the determination of the stel- 
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Figure 1. The stellar mass function of Virgo galaxies in the clus¬ 
ter core (the ‘pilot region’), excluding M87. The dotted (blue) line 
indicates detected objects; the solid (black) line includes a cor¬ 
rection for estimated incompleteness. The error contours indicate 
Poisson errors plus a fixed 10% error accounting for systematics in 
the stellar mass determinations. The dashed (red) line is an inte¬ 
grated Schechter function fit with M* = IO^^Mq, (p* = 5.5 and 
a = -1.35. 

lar mass function for the cluster core. In section [3] we 
describe our cluster simulations. In section |4] we discuss 
the subhalo catalogs, and present three different models 
for the infall mass distribution. In sectionjSjwe match the 
subhalo mass function to the observed stellar mass func¬ 
tion, and discuss uncertainties in the resulting SHMR. 
We summarize our results in section El In an appendix, 
we also review subhalo properties, and discuss the trends 
in infall redshift with radial separation or velocity off¬ 
set from the cluster center. Throughout t his paper we 
assume a WMAP-7 (|Komatsu et al.l [201111 ACDM cos¬ 
mology with rif, = o.odiQ, ric = 0 . 220 , Ha = 0 . 729 , 
h = 0.703, n = 0.966, and erg = 0.809. 

2. THE NGVS STELLAR MASS FUNCTION 

The NGVS is a multi-band, panoramic survey of the 
Virgo cluster using MegaGam at the Megaprime focus 
on CFHT. An introduction to the survey, data analysis 
pipelines, data products and science goals is given in Pa¬ 
per I. Here we summarize information about the survey 
relevant to this work. 

The survey covers 104 square degrees around the two 
main components of the Virgo cluster, the A and B sub¬ 
clusters, extending out to 1-1.5 times the projected virial 
radius, depending on the estimated mass, in four bands 
u*, g, i, and z in the MegaGam filter systen:13 (for simplic¬ 
ity we will refer to these as u,g,i, and z in the remainder 
of the paper). 

NGVS represents the first major update to our inven- 

^ Note we assume this baryon density to generate a realistic ini¬ 
tial power spectrum, although the subsequent numerical evolution 
treats all matter as collisionless. 

^ The survey was originally designed to include r-band coverage, 
although this is not yet complete outside the pilot region at the 
time of writing. 


tory of V irgo galaxies s i nce th e Virgo Gluster Gatalog 
(VCG) of lBinggeli et al.l l)1985h . The VCC covered most 
of the region surveyed by the NGVS, as well as part of the 
southern extension, with photographic plates in a single 
filter, reaching a depth of Hiim ~ 20 mag for galaxies 
or 22-23 mag for point sources, and a surface bright¬ 
ness limit of ~ 25 magarcsec“^. By comparison, NGVS 
reaches a limit oi g = 25.9 mag for point sources at S/N 
= 10, and a surface brightness limit of 29 magarcsec“^ in 
g. The exceptionally faint surface brightness limit is the 
result of a new strategy for removing scattered light and 
treating large-scale spatial variations in the background, 
Elixir-LSB, summarized in Paper I. The excellent see¬ 
ing conditions at the CFHT site (NGVS images have a 
FWHM in the i-band < 0.6") also allow us to resolve all 
but the most compact objects at the distance of Virgo, 
separating most ultra-compact dwarfs (UCDs) and even 
some of the more extended globular clusters from point 
sources. 

Given the timescale for obtaining, calibrating and an¬ 
alyzing survey data, the observations of an initial ‘pilot’ 
region were analyzed first and it is these data we con¬ 
sider in this paper. We will present a detailed analysis 
of the full cluster stellar mass function in a subsequent 
paper. The pilot region consists of four contiguous fields 
(4 square degrees total area) around M87, the central 
galaxy of the main component of Virgo. These fields 
have complete coverage in all five bands, both in long 
exposures and in the short exposures designed to sup¬ 
plement our information in the saturated areas of the 
brightest galaxies. 

The identification and characterization of Virgo cluster 
members will be discussed in detail in an upcoming pa¬ 
per of this series (Ferrarese et al. 2015, in preparation); 
here we summarize the most relevant points. The detec¬ 
tion of extended, low surface brightness Virgo galaxies 
in a field contaminated with a far larger number of fore¬ 
ground stars and background objects i s a challenge - con¬ 
venti onal codes, such as SExtractor (|Bertin fc Arnouti 
Il996f) . would always regard low surface brightness objects 
as belonging to compact and/or brighter contam inants. 
To ci rcumvent the problem, a ring median filter (iSeckerl 
II995I1 was first applied to each g-band image (the NGVS 
images with the highest signal-to-noise ratio), with ra¬ 
dius adjusted to suppress all unresolved sources (stars 
and globular clusters) as well as compact background 
galaxies. A specific optimization of SExtractor was then 
run on the medianed stacks, thus allowing the identifica¬ 
tion of low surface brightness sources that are potential 
Virgo candidates. The resulting object catalog does not 
of course discriminate between Virgo members and con¬ 
taminants; therefore the next challenge is to define mem¬ 
bership criteria. These are based on the location of each 
galaxy in a multi-parameter space defined by a combina¬ 
tion of galaxy structural parameters (specific ally size and 
surface brightness, measured using GalFit (|Peng et al.l 
[ 2003 ), in appropriately masked original images), pho¬ 
tometric redshifts (based on u,g,i,z and, when avail¬ 
able, r-band photometry), and an index measuring the 
strength of residual structures in images created by sub¬ 
tracting from each galaxy the best-htting GalEit model. 
The exact combination of axes in this space was selected 
to allow for maximum separation of known Virgo and 
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background sources. The former comprise spectroscopi¬ 
cally confirmed (mostly) VCC galaxies, while the latter 
are identified, using the same procedure described above, 
in four control fields located three virial radii away from 
M87 and presumed to be devoid of cluster members. 

The reliability of the procedure was tested in two in¬ 
dependent ways. First, approximately 40,000 artificial 
galaxies were injected in the frames spanning the full 
range of luminosity and structural parameters expected 
for genuine Virgo galaxies, and then processed as de¬ 
scribed above. With few exceptions (for instance galaxies 
that land in the immediate vicinity of bright saturated 
stars or near the cores of high surface brightness galax¬ 
ies), all galaxies with surface brightness high enough to 
be visible in the frames were indeed recovered. A detailed 
discussion of the completeness of the data and the biases 
in the recovered parameters will be included in Ferrarese 
et al. 2015 (in prep.). Second, three NGVS team mem¬ 
bers independently inspected the 4 square degree of the 
pilot project region and identified, by eye, all objects that 
appeared to be bona-fide Virgo members; again with few 
exceptions, all such galaxies were detected and correctly 
identified as Virgo members by the code. 

The photometry for each Virgo member was recovered 
using three s eparate techniques: by fitting a single Sersic 
(|Sersiclll963ll component as implemented by GalFit; by 
fitting a core-Sersic, Sersic or double Sersic (for galax¬ 
ies with a stellar nucleus) model to the 1-dimensional 
surface brightness profile derived by finding the ellipti¬ 
cal isophotes that best fit the galaxy’s surface brightness 
distribution; and by a non-parametric curve-of-growth 
analysis. There is generally excellent agreement between 
the magnitudes derived with the three different meth¬ 
ods. All magnitudes were dereddened using reddening 
estimate d at the locati o n of the center of each galaxy 
following iSchlegel et al.l (jl998l l. 

Tests with the artificial galaxies also allowed for a strict 
quantitative assessment of our detection and complete¬ 
ness limits and galaxy parameter recovery. For a given 
magnitude, g, we describe the completeness, f{g), as 
the fraction of input artificial galaxies recovered by the 
pipeline. Scatter and bias in the recovered magnitudes 
were investigated and turn out to have only minor impact 
on the final stellar mass function, whereas completeness 
corrections become significant below g = 19. Based on 
these tests, we fitted the recovered fraction as /(g) = 1 if 
g < 18.9 and f{g) = 0.511 - 0.2a;-k0.00444a;3 if g > 18.9, 
where x = g — 22. The true number of objects was then 
estimated as l//(g) of those detected. 

Stellar masses were obtained via Bayesian modeling 
of the ugriz spectral energy distributions (SEDs) of the 
Virgo members. The integrated magnitudes of these 
galaxies were measured via a curve-of-growth analysis 
of the multi-wavelength imaging from the NGVS. Errors 
were assigned to the photometry statistically by com¬ 
paring the curve-of-growth magnitudes to those from 
an independent GALFIT analysis of the same imag¬ 
ing. Our stellar population synthesis models span a 
multi-dimensional parameter space designed to mimic 
the wide variety of star formation histories (SFHs), chem¬ 
ical enrichment histories and dust contents of present- 
day galaxies. We empl oyed the base SSP models of 
iBruzual fc Gharlc^ (1200,81) and extinction was treated fol¬ 


lowing the two-component prescription of iGharlot &: Falll 
(|2000r) . A finite grid of 50,000 synthetic stellar pop¬ 
ulations were then genera t ed, a ssuming the priors de¬ 
scribed in Ida Gunha et al.l (120081 ). We fitted the SED of 
each galaxy using this grid to determine the marginal¬ 
ized posterior PDF for its stellar mass. The final value 
corresponds to the median value of this PDF, with an 
uncertainty equal to half the interval between the 16th 
and 84th percentiles. 

An important caveat regarding our modeling is the ef¬ 
fect of bursts of SF on the resultant M*/T. Our priors 
assume that bursts have a 50% probability of occurring 
at each timestep throughout the lifetime of our synthetic 
populations, and that half of the models have not ex¬ 
perienced a burst within the past 2 Gyr. Since young 
(i.e., bright) stars have lower M*/L, it would be fair to 
say that our fiducial models are biased to low M,/L. 
Reducing the contribution of bursts to the SFHs of our 
synthetic populations indicates that this bias could be as 
high as 0.2 dex. However, given that the priors of our 
fiducial model have been tailored to reproduce the spec¬ 
troscopic properties of SDSS galaxies, that the SFHs of 
real galaxies are unlikely to be smoothly varying, and 
that broadband colors cannot effectively constrain the 
role of bursts, we retain the M* values predicted by our 
fiducial model for this analysis. Further details on the 
SED modeling of NGVS galaxies will be presented in 
a forthcoming paper on the stellar populations of the 
Virgo cluster (Roediger et al. 2015, in preparation). 

Finally, we note that we have not included UCDs in 
the core stellar mass function, since it is unclear whether 
to treat these as independent star-forming systems, or a 
population o f unusually mass i ve clusters associat ed with 
galaxies (see iZhu et ahl 120141: iZhang et al.l 120151 : Liu et 
al. 2015, in preparation, for a detailed discussion and 
references). The several hundred UGDs in the core re¬ 
gion could potentially double our stellar mass function 
at the low-mass end, halving the halo masses derived for 
the smallest objects by abundance matching. Further 
work on the dynamics of these objects may clarify their 
cosmological status. 

Figure [T] shows the cumulative stellar mass function for 
the core region, excluding M87. The dotted (blue) line 
indicates detected objects; the solid (black) line includes 
the correction for incompleteness. The error contours 
indicate Poisson errors plus a fixed 10% error account¬ 
ing for systematics in the stellar mass determinations. 
(Using SFHs with fewer bursts could further shift the 
stellar mass function up to 0.2 dec to the right, produc¬ 
ing an increase in counts roughly equal to the l-cr error 
range shown here.) Overall, for the purpose of this work, 
our sample consists of 407 galaxies with reliable stellar 
masses in the core region, spanning five decades or more 
in stellar mass. Below M* ^ 5 x 10®Mq, the cumula¬ 
tive mass function is roughly a power law. An integrated 
Schechter function (ISF) 

pOO 

N{>M)= (j)*x°^exp{-M/M*)dM/M\ 

Jm 

with M* = IO^^Mq, (j>* = 5.5 and a = —1.35 provides 
an excellent description of the stellar mass function (ex¬ 
cluding M87) everywhere except at the smallest masses. 
Below M* = 3 X IO^Mq, the mass function shows a pos- 
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sible change in slope, but the reality of this feature is 
unclear given systematic uncertainties in the complete¬ 
ness corrections at the faintest magnitudes. Thus in what 
follows we will adopt the ISF form with the parameters 
given above to represent the stellar mass function. 

The ISF has degeneracies in the model parameters, 
and we did not find that automated determination of 
all three parameters simultaneously worked particularly 
well. Instead, we determined experimentally that for 
fixed M* = IO^^Mq, a plausible range of slopes is 
OL = (—1.3)-(—1.42), and that the normalization can be 
in the range (jf = 3-8, depending on the value of a. 
Allowing for SFHs with fewer bursts and systematically 
higher mass-to-light ratios broadens this range slightly, 
permitting a flatter fit with a (—1.28) and <jf = 12. 
In Figure |9] below, we consider several extreme choices 
of ISF parameter values that still provide a reasonable 
approximation to the data. 


3. NUMERICAL SIMULATIONS 

Abundance matching requires an estimate of the num¬ 
ber of dark matter structures in a region, either from 
analytic theory or from numerical simulations of struc¬ 
ture formation. In our case, the region in question 
corresponds to a line of sight through the core of a 
cluster, so most of the galaxies we detect should oc¬ 
cupy subhalos within the larger cluster halo. While 


dance do exist fe. 


Tavlor & Babull 12004 I2005allbl: 

van den Bosch et al. 

120051: IGiocoli et al.ll2008l: IGan et alJ 

2010HYang et al.ll2011l). their accuracy is unclear in the 


densest environments. Thus, we will use numerical N- 
body simulations of halo formation to estimate the num¬ 
ber of subhalos seen at small projected separation from 
the cluster core, and to determine their characteristic 
properties, e.g., their infall redshift, or their original mass 
at the time of infall. Our simulations include dark mat¬ 
ter only, since the baryonic physics affecting the cluster 
mass distribution is complex and mo del-dependent. We 
discuss several specific baryonic processes which may af¬ 
fect our results in section U . 

In order to leverage the full power of the NGVS obser¬ 
vations, we need simulations that resolve substructure 
down to very small masses. To accomplish this, we use 
the technique of resimulation, in which a high-resolution 
region of interest is embedded within a lower-resolution 
simulation. Multi-resolution initial conditions were gen - 
erated using the Grafic2 package l|Bertschingerl 120011). 
and e volved using the V-body code Gadget2 (|SDringell 
120051) . The initial, top-level simulation had 256^ particles 
in a cubic volume 140 h~^ Mpc on a side, large enough 
to contain many Virgo-mass clusters. With the chosen 
volume and WMAP-7 cosmology, this corresponds to a 
mass of 1.23 x 10^° h~^MQ per particle. The softening 
length used was 0.02 times the mean inter-particle sepa¬ 
ration, or e = 10.94 h~^ kpc in comoving units. In the 
top-level simulation, CDM halos were detected using the 
University of Washington friends-of-friends code FOP0. 
Gandidate halos for resimulation were selected from the 
mass range 2-5 xIO^^Mq, comparable to the estimated 
mass of the Virgo cluster (Paper I). In order to mini¬ 
mize interference from other large halos, any halo with 


® http://www-hpcc.astro.washington.edu/tools 


a neighboring halo of greater than i/s its mass within 
3 h~^ Mpc was excluded from consideration. From the 
halos meeting these criteria, 10 were chosen at random 
for resimulation. These 10 halos represent a wide range 
of formation epochs and mass accretion histories. 

Each of the selected halos was resimulated individually. 
For each halo, we determined the smallest rectangular 
volume in the initial conditions that contained all the 
particles in the final halo at z = 0. This rectangular 
region was then extended to twice the linear size in each 
dimension, and all particles within this larger volume 
were replaced with higher-resolution initial conditions, 
using Grafic2. By choosing a larger volume around the 
cluster, we ensure that all particles associated with the 
final halo, out to roughly twice the virial radius, will 
be high-resolution particles, and that none of the more 
massive particles in the larger volume will intrude into 
the halo proper. 

The resimulations were performed with a factor of 
1000 increase in mass resolution, the mass of the high- 
resolution particles being 1.23 x lO"^ For the 

high-resolution particles, a softening length of 0.014 
times the mean (high-resolution) inter-particle separa¬ 
tion was used, Chigh = 0.7656 h~^ kpc in comoving units. 
This choice corresponds to the optimal softening length 
for resolving substr ucture in halos of t he selected mass 
range, as defined in iPower et al.l (|2003f) . The softening 
length used for the low-resolution particles was the same 
as in the top-level simulation, eiow = 10.9375 h~^ kpc. To 
record the assembly of each cluster we output 121 snap¬ 
shots, equally spaced in logo betw een z = 0 to 9. We 
fit Navarro-Frenk-White (NFW - iNavarro et al.l [19971) 
profiles to the main cluster to determine its total mass, 
virial radius and concentration; neighboring halos and 
cluster subhalos were also detected in each snapshot, as 
described below. The structural properties of the resimu¬ 
lated halos are listed in Table [T] (Masses and virial radii 
assume the spherical collapse definition of the virial over¬ 
density, rather than a fixed density contrast such as 200 
times the critical density pc or the matter density pm-) 
They contain from 8.8 to 24.4 million particles, equiv¬ 
alent to a mass of 1.54-4.28 xIO^^Mq, and have sev¬ 
eral thousand resolved subhalos and sub-subhalos within 
their virial volume. 

The clusters in our final sample vary significantly in 
mass, formation history and concentration, as indicated 
in Table [H We could restrict the sample further, choos¬ 
ing a particular range of concentrations and/or formation 
histories based on more detailed dynamical or structural 
modeling of Virgo. Instead, we will retain the entire sam¬ 
ple and use the cluster-to-cluster scatter as an indication 
of the systematic uncertainties in the subhalo mass func¬ 
tion introduced by variations in the global structure and 
history of the cluster. We will normalize our results to a 
common cluster mass, however, since subhalo mass func¬ 
tions scale roughly with the mass of the parent halo. We 
take the mass of Virgo to be M 2 oo,c = 4.2 x lO^^Mr:^ and 
its con centration to be c = 2.51, following iMcLaughlinI 
(119991) . Assuming an NFW profile with the correspond¬ 
ing scale radius, this corresponds to a spherical col¬ 
lapse mass Mvir = 5.76 x IO^^Mq and concentration 
Cvir — 3.5. These are the fiducial values adopted by 
the NGVS survey for the main component of Virgo, 
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Halo 


Afpart (10«) 

.^sub 

rvir (kpc) 

np (kms ^) 

c 

Z50 

3 

3.67 

20.9 

5328 

1868 

954 

4.1 

0.42 

6 

3.66 

20.9 

4983 

1867 

986 

5.4 

0.69 

7 

2.04 

11.7 

2366 

1536 

907 

9.8 

1.30 

15 

2.03 

11.6 

2546 

1534 

827 

6.1 

0.99 

20 

2.39 

13.7 

3437 

1619 

843 

4.9 

0.67 

24 

2.69 

15.4 

3656 

1685 

883 

5.1 

0.52 

27 

1.54 

8.8 

2154 

1400 

841 

10.6 

0.39 

29 

4.28 

24.4 

5786 

1967 

1071 

6.5 

0.45 

35 

1.86 

10.6 

2009 

1489 

819 

6.9 

0.66 

40 

2.54 

14.5 

3756 

1653 

862 

4.9 

0.64 


Table 1 

Properties of the resimulated halos: total halo mass number of particles A^part, number of subhalos A^sub? virial radius rvir» peak 
circular velocity Up, NFW concentration parameter c, half-mass accretion redshift 250 


component A. Recent work bv lUrban et alJ (|20lH) . us¬ 
ing X-ray spectroscopy and scaling relations, found a 
mass three times smaller for component A. It is also un¬ 
clear whether component B, the other large component 
of Virgo, has merged with component A or should be 
considered a distinct halo. Component B has a mass of 
Af200,c 1 X 10 ^"^Mq (Paper I), so this would increase 

the mass of the cluster by 20% . Thus the mass of the 
cluster could conceivably lie in the range 0.33-1.2 times 
our fiducial value. We will discuss the implications of a 
different mass for Virgo in section [5j 

4. SUBHALO CATALOGS 

From our resimulated halos, we then created subhalo 
catalogs to match to the NGVS observations. A rich va- 
riety of halo-finding techniq ues exist in the literature (cf. 
iKnollmann &: Knebell200^ for a review), but only a few 
are suited to finding subhalos within larger bound struc¬ 
tures. To find subhalos in the z = 0 snapshot, we used 
the Amiga Halo Finder (AHF - IKnollmann fc Knebel 
IM)^ . For the preceding snapshots, which were used 
only to identify the infall redshift and infall mass for 
each subhalo, we detected halos using the simpler and 
faster FOF code. We describe each of these steps below. 

4.1. The Subhalo Mass Function at z = 0 

In the z = 0 snapshot, we need complete information 
about subhalos, including their position, velocity, final 
mass, and internal substructure. To derive the subhalo 
catalog we used the public halo finder AHF. AHF is a 
multi-scale, hierarchical group finder designed to run in 
parallel on large simulations, using the message-passing 
protocol (MPI). It constructs an initial representation of 
the density distribution sampled on a coarse grid, iden¬ 
tifies isolated regions that exceed some density thresh¬ 
old above the mean, and refines the grid iteratively in 
these regions. Structures are then identified as density 
peaks within this hierarchy, either at a base level in the 
field (halos) or as substructure (subhalos) within larger 
structures. As part of the process of associating par¬ 
ticles with peaks, unbound particles are iteratively re¬ 
moved from each object, to ensure that it represents a 
bound physical structure. To identify halos at z = 0, we 
use AHF’s built-in calculation for the (spherical collapse) 
virial overdensity, rather than a fixed overdensity of 200 
Pc or Pin. Hereafter, we will refer to the basic subhalo 
catalog produced by AHF as ‘model O’. 

http://popia.ft.uam.es/AHF 
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Figure 2. The mean subhalo (self-bound) mass function at red¬ 
shift 2 = 0, determined using the group finder AHF. Dotted lines 
show the 1-cr simulation-to-simulation scatter. The lower axis 
shows subhalo mass at redshift zero, Mq, relative to the virial 
(spherical collapse) mass of the main halo Mmain- The upper axis 
shows subhalo mass normalized to our fiducial (spherical collapse) 
mass for Virgo, Mmain = Afvir = 5.76 x 10 ^'^Mq. 

Figure [2] shows the cumulative subhalo (self-bound) 
mass function obtained from AHF at z = 0, using the 
AHF parameter value NperRefCell = 4, and a minimum 
group size of 20 particles. The thick solid line is the 
average for the sample of 10 clusters, while the dot¬ 
ted lines show the 1 -ct cluster-to-cluster scatter. The 
bottom axis shows subhalo mass at redshift zero. Mg, 
relative to the mass of the main halo, while the top 
axis shows Mg^eqv, the subhalo mass at redshift zero 
rescaled to our fiducial spherical collapse mass for Virgo, 
Mvir = 5.76 X lO ^^Mf;^. Consistent with many pre¬ 


vious results (e. 

e.. IDiemand et al.l 

20071 ISorinEel et al. 

2008; 

Gao et al. 

2011 

: iKlvoin et al. 

120111: IGontini et al. 

2012 ; 

Gao et al. 

2012 

), we find an approximate power- 


law with a logarithmic slope din TV (> M)/d\nM ~ —0.9 
over almost four decades in mass, extending to a max¬ 
imum mass ratio of Mg/Mmain “^2-5% . The slope be¬ 
comes slightly shallower around A//Mmain ~ 2-3x10“® 
(Meqv ^ 10^Mq), where resolution effects become im¬ 
portant. 
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4.2. Tracing Subhalo Histories 

The higher-redshift snapshots were used mainly to con¬ 
struct histories for each AHF subhalo detected at 2 ; = 0, 
to determine when they last merged with a larger sys¬ 
tem. To construct the required merger tree^ we used the 
simpler and faster group-finding code FOFQ, modified to 
work with multi-scale Gadget2 snapshots. The standard 
linking length of 0.2 times the mean inter-particle sep¬ 
aration and a minimum group size of 20 particles were 
used to define the FOF groups. 

We note there is a slight difference in our mass defi¬ 
nitions in different time steps, since the AHF results at 
z = 0 are spherical overdensity (SO) masses , whereas the 
result s at higher redshift are FOF masses. I Tinker et al.l 
(j2008H compared the two types of mass estimates, and 
found that they agree to with 5-10% on average, although 
SO masses can be up to 30% smaller in cases where FOF 
artificially links nearby structures. AHF avoids these ex¬ 
tended structures by using an additional u nbinding step 
in its m ass estimates. We note however that iTinker et al.1 
(j2008[l uses a SO method with an overdensity threshold 
of 200pm, which produces masses up to 10% higher than 
the AHF threshold at low redshift. As a result, our es¬ 
timates of Minf/Mmain may be 5-10% larger than values 
calculated entirely using AHF, for systems merging at 
z ^ 0.5. Given that we are interested in broad trends 
over many orders of magnitude in mass, and that rela¬ 
tively few subhalos merge at redshifts this low, we will 
ignore this offset. 

To determine the merger history of each AHF subhalo, 
we proceed through the following steps: 

1. First, we record all the subhalo’s particles in the 
final snapshot, and search the preceding snapshot 
for structures (halos or subhalos) containing these 
same particles. We define the ancestor of the sub¬ 
halo as the structure containing the largest fraction 
of its particles in the preceding snapshot. 

2. Iterating backwards through the snapshots from 
ancestor to ancestor, we determine an evolutionary 
history for each subhalo. We distinguish parts of 
the history where the ancestors are ‘merged’ sub¬ 
halos in a larger halo from parts where they are 
‘independent’ halos in the field. 

3. We define the maximum mass, as the largest 

mass a subhalo’s ancestor ever had as an indepen¬ 
dent halo. The infall redshift Zinf (or infall time 
tinf) is taken to be the last redshift (time) at which 
the ancestor existed as an independent halo before 
becoming part of a larger system, and the infall 
mass Mini is defined as the mass of the ancestor 
in the final snapshot before it became part of the 
larger system. 

4. At each step, the ‘main halo’ is defined to be the 
ancestor of the main halo in the subsequent step, 
that is it is always the halo that contains the largest 
fraction of the particles of the main halo of the sub¬ 
sequent step. By tracing the main halo back step- 
by-step, we can define a ‘main trunk’ of the merger 

® See footnote 3. 


tree, as opposed to the smaller ‘side branches’ that 
merge with it. 

This analysis assumes that merging is a well-defined, 
one-way event. In fact, FOF often links together, in the 
early phases of merging, structures which then separate 
back out into distinct halos in later st eps. This ‘false 
merge r’ problem is discussed in detail in iFakhouri fc 
(|2008H . To reduce the influence of false mergers on our 
final statistics, we proceed forward in time through each 
subhalo’s history and save the current state whenever its 
mass increases by more than 1/3 between snapshots, indi¬ 
cating that a major merger has taken place. If the mass 
subsequently decreases by more than 1/4, this almost cer¬ 
tainly represents the end of a false merger, so we revert to 
the last saved state and exclude from the calculation of 
the maximum mass all ancestors between that snapshot 
and the snapshot under consideration. In this way, tem¬ 
porary and artificial jumps in mass from false mergers 
do not incorrectly reset the maximum mass. 

It is not always possible to find Minf and Mmax, espe¬ 
cially for the smallest halos. It can be that a subhalo’s 
ancestor was never independent, for instance, given our 
time resolution and starting redshift. This is not a signifi¬ 
cant problem, however, affecting less than 1% of subhalos 
overall, and almost none of those with ~ 200 particles or 
more. 

4.3. Correcting for Completeness and Numerical Effects 

Three sorts of problems can arise when creating sub¬ 
halo catalogs, that cause them to be incomplete. First, 
the group-finder may not correctly identify substructure 
in a given simulation output, or may not associate with 
a subhalo all the particles that are actually bound to 
it. This is especially likely for low-mass sub halos and/or 
subh alos in the dense core of the host halo (|Knebe et al.l 
Second, due to the discrete, iV-body nature of 
the simulations, the dissipation of substruc ture is accel¬ 
erated by numerical relaxa tion effects (e.g.,lMoore et al.l 
II996I iDiemand et al.ll2004 see IMoor^ 120001 for a histor¬ 
ical review of this problem). Thus the subhalos in the 
final snapshot will generally be less massive than they 
would have been if simulated at higher resolution, and 
some subhalos that should have survived to the present 
day will end up being completely disrupted. This prob¬ 
lem affects the least massive and oldest subhalos most 
strongly. Third, there is an important difference be¬ 
tween the simulations and the physical process being 
modeled, namely the presence of baryons. Generically, 
baryons should cool and condense into the cores of sub¬ 
halos, making them more resistant to mass loss and 
disruption. On the other hand, in some cases feed¬ 
back from star formation or nuclear activity may eject 
mass from the center of the system, helping to dis¬ 
rupt the subhalo. Overall, we expect baryonic simula¬ 
tions to retain more substructure on galactic scales than 
the collisionless simulations co nsidered here, but the de¬ 
tails will be model-dependent (IRomano-Dfaz et al.ll2009L 
l20ICtl^hewtschenko fc Macci6ll20lTlb 

While an alternative choice of gr oup-finder might re - 
cover some of the missing structure (|Knebe et al.ll20l3l , 
numerical relaxation and the lack of baryons in the 
simulation are harder problems to address. Including 
baryons in the simulation brings with it the attendant 
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problem of simulating galaxy formation correctly; and 
although overmerging can be decreased by increasing 
the resolution of the simulation, convergence is expected 
to be relatively slow , especially in the central regions 
(jDiemand et al.l[2004fl . Instead, we will consider alterna¬ 
tive approaches to recovering structure, based on subhalo 
merger histories. 

4.4. Three Estimates of the Infall Mass Function 

4.4.1. Model 1 - An Upper Limit on the Mass Function 

To overcome the issue of incompleteness with the sub¬ 
halo catalogs, we can try a different approach. Every 
subhalo in the main halo was at one time an indepen¬ 
dent halo in the field. Thus, every subhalo will corre¬ 
spond to a merger event in the merger tree. We can 
therefore use these mergers as a way to identify subhalos. 
Similar ‘historical’ approaches to finding substructure 
have been implem ented previously with several group - 
finding codes le.g.. IGill et al.l[2004 iTormen et al.l 20041 . 
In the most conservative limit, we assume that every 
halo merging at any point in the merger tree, whether 
with the main trunk or with a side branch, survives to 
the present day as a self-bound substruct ure. This ap¬ 
proach is similar to those u s ed by , e.g., lTormen_^t_aL 
( 2004 ). Ivan den Bosch et'all (|2005[) . and IGiocoli et al. 

( 20081 ) . except that we include side branches, account¬ 
ing for higher-order substru cture (sub-subha l os), a s in 
the semi-ana l ytic m odels of iTavlor fc Babull (120041) or 
IGiocoli et al.l (j20inD . In order to obtain z = 0 subhalo 
properties, such as position and velocity, we further as¬ 
sume that the subhalo is traced by the most bound par¬ 
ticle of its last independent ancestor. It is, of course, 
impossible to obtain z = 0 subhalo masses with this 
method, but we can measure the mass at infall Mjnf for 
each subhalo. Assuming that every merger has a cor¬ 
responding subhalo that survives down to z = 0, we are 
certain to capture all actual subhalos down to some mass 
limit, so this approach (‘model 1’ hereafter) provides a 
conservative upper limit on the SIMF at z = 0. Of 
course, model 1 also includes many systems that would 
not survive as distinct subhalos , since some will merge 
together physically or be disrupted, either before or af¬ 
ter they enter the main halo. If, however, we compare 
model 1 to the previous results from AHF (‘model O’, de¬ 
fined above), which detects almost everything that does 
survive despite the factors that artificially accelerate the 
disruption of substructure in our simulations, we can 
constrain the SIMF both from above and from below. 

4.4.2. Model 2 - An Intermediate Model 

Model 1 makes several unrealistic assumptions. First, 
it assumes that every system that merges with the main 
trunk of the merger tree represents a single, well-defined 
and dynamically relaxed halo. In the simulation, merging 
groups may have internal substructure or separate com¬ 
ponents that survive the merger as a distinct subhalos. 
Physically, these systems might correspond to individual 
galaxies in a loose group that are dissociated from one- 
another when they merge with the cluster. We can cor¬ 
rect for this effect in part by identifying the least relaxed 
systems as they merge, and splitting them into their sub¬ 
components. To do this we use the relaxation parameter 
Xoa = |aicOM — ^MBl/j’vir, where fcOM is the position of 


[Mel 
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Figure 3. The mean subhalo infall mass mass function (SIMF), 
using the three different models described in the text. The dashed 
(red) curve includes only the subhalos found by AHF (model 0). 
The dotted (blue) curve includes every object that ever merged into 
the merger tree (model 1), while the solid (black) curve attempts to 
correct for loosely bound substructure and substructure disruption 
(model 2). The upper axis shows subhalo mass rescaled to the 
fiducial mass of Virgo, as in Figure (2] 

the group center of mass and xmr is the position of the 
most bound particle (jMaccib et al.ll200^ . If XoS > 0.25, 
we consider the merging subhalo sufficiently unrelaxed 
that it will fall in as one or more distinct components, 
and/or be tidally disrupted into these components by the 
present day. (a;off = 0.25 is approximately the relaxation 
parameter expected for a 3:1 merger where the secondary 
has just crossed the virial radius of the primary.) 

Second, we should also correct for the opposite limit, 
in which a sub-component of a subhalo remains bound 
to its parent subhalo after infall, and has merged with 
its parent completely by the present day. We find this 
is mainly an issue for massive subcomponents of sys¬ 
tems merging at z < 1. From dynamical friction ar¬ 
guments we expect merger or disruption rates to be 
fastest for sub-subhalos with a large mass relative to their 
parent. In particula r, seini-analytic calculations (e.g., 
ITavlor fc Babull I2004D suggest that the critical mass ra¬ 
tio for rapid disruption is around Msub/Afmain ~ 0.1. In 
practice, we consider that systems with Msub/Mmain > 
0.07 will have merged if tjnf is less than half the age of the 
universe at z = 0. The value 0.07 is chosen to match the 
model 0 results at the high-mass end of the subhalo mass 
function, where we expect the AHF substructure catalog 
to be complete and to provide good mass estimates. 

Although these two corrections are approximate, they 
prune the model 1 results down to a more realistic set 
of subhalos. Applying the two corrections to the merger 
tree, we define a ‘model 2’ for the infall mass function 
that resembles model 0 (the raw AHF results) at the 
high-mass end, and is intermediate between models 0 
and 1 at the low-mass end. 

4.4.3. The Three Models Compared 

































NGVS. IX. Abundance Matching in the Virgo Cluster 


9 


The three substructure models are compared in Fig¬ 
ures H II and [5] Figure |3] compares the mean infall 
mass functions predicted by the three models. Roughly 
speaking, the predicted form of the SIMF is a power-law 
with an exponential cutoff at Minf/Mmain 0.1. The 
power-law slope a = d\nN/dlnM is ~ —0.8 for model 
0, ~ —0.96 for model 1, and ~ —0.87 for model 2. These 
values are consistent with previous numerical and semi- 
analytic estimat es of the SIMF, which is expected to be 
fairly universal (|van den Bosch et al.ll2005l : iCiocoli et al.l 
12001^ . Methods that only count ‘first-order’ (i.e., di¬ 
rect) branchings off the main trunk of the merger tree 
find slopes of ^ ~0 - 8 or s hallower fe.g.. iTavlor fc Babull 
12004 ICiocoli et al.ll2008ll . while methods that account 
for higher-orde r branchings (or sub-subhalos) find slopes 
of ^ —0.9 te.g lT avlor fc Babull[2004 ICiocoli et al.ll2010l : 

iJiang fc van den Boschll2014ll . Although it was not cho¬ 
sen to match these estimates, model 2 is quite close to 
the slope and norm alization derived in the most recent 
semi-analytic work (| Jiang fc van den Boschl[20l4l . 

Overall, the range in uncertainty between the three 
models is roughly a factor of 2-3 except at the small¬ 
est masses, where numerical resolution renders the raw 
AHF results (model 0 - lowest curve) increasingly in¬ 
complete. Model 2 predicts results similar to model 0 at 
large masses, but finds a factor of 2 or more subhalos at 
small masses. Figure |4] shows subhalo abundance as a 
function of infall mass and infall redshift. We see that 
both models 1 and 2 (right and center panels) include a 
large number of early and/or low-mass subhalos, relative 
to model 0. On the other hand, all three models find 
similar numbers of massive, recently merged halos. 

Since older, more stripped subhalos will generally lie 
on more tightly bound orbits, we also expect the three 
models to predict different radial and velocity distribu¬ 
tions. Figure [5] shows subhalo abundance as a function 
of radial offset from the cluster center (left panel), and 
line of sight velocity offset (right panel), normalized to 
the virial radius and the circular velocity at the virial 
radius, respectively. There is a dramatic difference in 
clustering between model 0 (the raw AHF results) and 
models 1 and 2, but relatively little difference in cluster¬ 
ing between models 1 and 2. 

Thus, in summary, averaged over the whole virial vol¬ 
ume model 1 predicts roughly twice the abundance of 
substructure at a fixed infall mass, while model 2 pre¬ 
dicts numbers closer to model 0, but models 1 and 2 
both have many more old, low-mass, and/or centrally lo¬ 
cated subhalos. When we restrict the subhalo sample to 
the central region of the cluster, the result is that the 
three models predict significantly different subhalo mass 
functions in this region. 

To compare the simulation results to the stellar mass 
function in the pilot region, we define an analogous re¬ 
gion in each simulation: a beam with a square cross- 
section 0.252 Tvir on a side, passing through the center 
of the cluster along the line of sight. For our fiducial 
Virgo mass, this corresponds to the part of the cluster 
covered by the pilot region , as discussed in section [5] 
below. Figure |6] shows the SIMF for halos in this cen¬ 
tral region. Within the central region the mass functions 
for models 1 and 2 now differ by a factor of more than 
5, highlighting the considerable uncertainties in subhalo 
abundance matching that are associated with the details 


of subhalo evolution. We will discuss the effect of these 
differences on abundance matching below; we also dis¬ 
cuss the spatial and orbital properties of the subhalos 
further in appendix lAl 

5. THE STELLAR-TO-HALO MASS RATIO 
5.1. Basic Results in the Cluster Core 

To determine the subhalo SHMR, we use the tech¬ 
nique of abundance matching, assuming the most mas¬ 
sive satellite galaxy is hosted by the most massive sub¬ 
halo, the next most massive by the next most massive 
subhalo, and so on. There are three definitions of the 
subhalo mass that we can use to construct this relation¬ 
ship: Mq, the mass of the subhalo at z = 0; Minf, the 
final mass the subhalo had as an independent halo, imme¬ 
diately prior to merging into a larger system; and Mmax, 
the largest mass the subhalo ever had as an indepen¬ 
dent halo. In practice, we find the results are almost 
identical for the latter two choices. Using Mq produces 
very different results (systematically lower halo masses 
for a given stellar mass), but this choice seems unphysical 
since tidal stripping may have reduced the dark-matter 
mass of a subhalo significantly without affecting its inner 
stellar component. Thus we consider Minf as the most 
logical choice for subhalo abundance matching. In par¬ 
ticular, if the SHMR were independent of redshift, then 
we would expect it to be comparable in the clusters and 
in the field, when expressed in terms of Minf. (There is 
in fact evidence that the field SHMR varies slightly with 
redshift, as discussed below.) 

The fiducial model of the Virgo cluster adopted by 
NCVS has a total mass M 200 .C = 4.2 x IO^^Mq and a 
concentration c = 2.51 with respect to the outer radius 
^ 200 ,c = 1-55 Mpc, or 5.38*^ at the fiducial distance of 
16.5 Mpc. Assuming a NFW profile, the corresponding 
spherical collapse mass is Mvir = 5.76 x IO^^Mq and the 
spherical collapse virial radius is Tvir = 2.16 Mpc, or 7.5° 
at the distance of Virgo, giving a spherical collapse con¬ 
centration Cvir 3.5. The pilot region corresponds to a 
square patch 2° on a side (4 MegaCam pointings), cen¬ 
tered on M87, the center of component A. Taking into 
account chip boundaries and edge effects, the effective 
area is 1.9° on a side, or 0.252 rvir on a side for our hdu- 
cial mass and distance. Thus, in each of three orthogonal 
projections of the resimulated halos, we select only the 
subhalos in a projected region of this size for comparison 
with the core stellar mass function, scaling the size of the 
region to the virial radius of each individual cluster. As 
explained in section O we also rescale the mass of each 
cluster halo to our fiducial mass for Virgo, and adjust all 
the subhalo masses accordingly. The 30 sets of subhalos 
(that is three projections of each of 10 clusters) thus de- 
Hned for each of the three substructure models comprise 
our simulation data for all of the following analysis, un¬ 
less otherwise stated. We match these to the previously 
determined observed stellar mass function (see section 
El), excluding the central galaxy M87, since it should be 
matched to the main halo, rather than any individual 
subhalo within it. 

The average SIMF for the central region, for the three 
different substructure models, was shown in Figure [HI 
Using the results for the intermediate model 2, which 
represents our best guess at the true mass function, we 
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Figure 4. Subhalo abundance as a function of infall mass and infall redshift Zinf. The left and right-most panels show the results 

for model 0 (AHF) and model 1, respectively. The central panel shows the results for the intermediate model 2. 



Figure 5. The normalized distributions of projected radial separation from the cluster center (left panel) and velocity offset (right panel) 
for the three substructure models; colors and line-styles are as in Figure!^ 


match the SIME for each of the 30 subhalo catalogs to 
the stellar mass function of the cluster core shown pre¬ 
viously in Eigure[TJ The results are shown in Figure 0 as 
total (infall) halo mass versus stellar mass (top panel), 
or halo-to-stellar-mass ratio versus stellar mass (bottom 
panel) . The thick solid (red) line shows the mean rela¬ 
tionship, while the thin dotted lines show the l-cr scatter 
amongst the 30 subhalo catalogs. The smooth curve at 
high stellar mass show s the results from the analysis of 
iLeauthaud et al.l (|2012l - L12 hereafter), which combines 
constraints from galaxy-galaxy leasing, galaxy cluster¬ 
ing, and the stellar mass function, for comparisorQ- 

We note that the results of Leauthaud et al. are based 
on isolated field galaxies, whereas we are considering a 
population of galaxies now for the most part located deep 
in the core of a galaxy cluster. On the other hand, we 
are using the subhalo infall mass Mi^i, rather than the 
present-day mass Mq, for comparison. If the present-day 

® Leauthaud et al. define halo mass using a density contrast of 
200 times the matter density pm, so their halo masses will be a tew 
percent larger than ours for identical halos at z = 0.88. 


cluster members were typical field galaxies at the mo¬ 
ment they fell in to the cluster, and if their stellar mass 
remained unchanged subsequently, we would expect their 
stellar-to-infall-mass ratio to match the SHMR for field 
galaxies at their infall redshift. Galaxies that fall into a 
cluster may not be completely typical, of course, as the 
conventional picture of structure formation predicts that 
they will always have resided in regions denser then the 
mean, even before they merged with the cluster. Com¬ 
paring our stellar-to-infall-mass ratio to the field SHMR 
allows us to explore this possible environmental depen¬ 
dence. The mean redshift for the L12 field sample 
shown here is (z) = 0.88, which should be slightly lower 
than the mean infall redshift of the subhalos, (zinf) = 
1.5 for model 2. We also note that L12 only measured 
the SHMR down to log(M^) ^ 11.5 at this redshift; the 
dotted line shows an extrapolation of their power-law fit 
to lower masses, to aid comparison. 

We see two main features in our results. First, at in¬ 
termediate stellar masses our results agree roughly with 
those of Leauthaud et ah, both in slope and in normal- 
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Figure 6. The average SIMF of the central region for the three 
substructure models. As in Figure 3, the dashed (red) curve is 
for model 0, the dotted (blue) curve is for model 1, and the solid 
(black) curve is for model 2. 


ization. There may be a slight 20%) offset in nor¬ 
malization, but it is only marginally significant given the 
cluster-to-cluster scatter and the difference in redshift of 
the two samples, and is also well within the range of 
possible bias in our stellar mass-to-light ratios discussed 
previously in section O Physical processes such as star 
formation or stripping of stellar material after infall could 
also affect stellar masses of cluster galaxies relative to 
the field; we would need better statistical evidence to 
establish whether this is the case. 

Combining our results at low stellar masses with those 
of L12 at high stellar masses, it appears the SHMR can 
be fit by a single power-law Mh from masses 

of to 10 ^°Mq (although at that 

point the slope of t he relationship starts to change). 
iBehroozi et al.l (|201Clll have proposed a more precise func¬ 
tional form for the relation between halo and stellar mass, 
also used in LI2: 


logio(M;,) = 


logio(Mi) -b/3 logio 




Over the range log(M*/M 0 ) = 5.0 - 10.5, our re¬ 
sults are well-described by this form with parameters 
[log(Mi/M 0 ), log(M*,o/M 0 ), /3, (5, and 7] = [12.45, 
10.35, 0.39, 0.4, 1.0]. Civen the uncertainties at the high- 
mass end, however, our data only really constrain the 
faint-end slope /3 and the normalization (set by Mi and 
M*.o)- Overall, measuring the abundance of faint dwarf 
galaxies in Virgo allows us to extend previous measure¬ 
ments of the SHMR several orders of magnitude in stel¬ 
lar mass, down to the faintest galaxies known outside the 
Local Group. 

The second feature is a possible offset between our re¬ 
sults and those of LI2 at large stellar masses. For the 



Figure 7. Halo mass Mh{top panel), and the ratio of halo to 
stellar mass Mfi/Mt, (lower panel), as a function of stellar mass 
M*. The thin lines show the l-tr scatter of the simulations. The 
smooth blue curves at large stellar mass show the results from 
ILeauthaud et ahl Il2012l ~l at 2 = 0.88 for comparison (solid, with an 
extrapolation of their fit to lower masses dotted). 


few largest objects in the core of Virgo, our results sug¬ 
gest infall halo masses ^^2-3 times larger than the results 
of L12 at (z) = 0.88. This could be because the largest 
systems merged at higher redshift, when the SHMR was 
lower in this mass range (see Figure [9] below), or it could 
indicate star formation has been suppressed in these sys¬ 
tems, relative to the field. Bias in our SHFs could also 
shift our stellar mass scale up (i.e., to the right in Figure 
|7|) by as much as to 0.2 dex, as discussed in section O 
The statistics are poor at this end of the relationship, 
however, due to the small number of bright galaxies in 
the pilot region. Results from the full NCVS survey will 
allow us to explore this trend over the full range in stellar 
mass. 


5.2. Comparison to Recent Work 
IBehroozi et al.l (|2013l . - B13 hereafter) recently mod- 
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eled average galaxy star formation histories and esti¬ 
mated the SHMR for field galaxies over a broad range 
of mass and redshift. Assuming our results within the 
core of Virgo reflect the field SHMR at the infall red- 
shift, they complement those of B13, covering most of 
the same redshift range, but extending down one to two 
orders of magnitude in halo mass at intermediate red¬ 
shift. 

Eigure [ 8 ] compares our derived SHMR to the results 
of B13 and L12. The thin colored lines show the mod¬ 
els of Behroozi et al. for z = 0.1, 1, 2, 3, 4 and 5 
(black, blue, cyan, green, yellow, and orange respec¬ 
tively), while the dashed lines show the results of Leau- 
thaud et al. for z = 0.37 and 0.88 (black and blue lines 
respectively). The thick black line shows our derived 
SHMR, and the dotted lines show the l-cr cluster-to- 
cluster scatter. Once again, our results lie slightly 
above the field SHMR at z=0-l, but are in fairly good 
agreement with the field estimates at redshift z ^ 2 , for 
the mass range 10^Mq- There is a possible 

discrepancy for M* > 3 x 10 ^°Mq, where we predict a 
factor of ^ 2 deficit of stellar mass associated with the 
most massive halos that merged into Virgo. Some or all 
of this offset could be related to the possible bias in our 
mass-to-light ratios discussed in section [51 

These preliminary results seem promising, but are sub¬ 
ject to a number of uncertainties. Figure [9] illustrates a 
few of the main uncertainties in the modeling. The top 
left panel compares our mean SHMR as a function of halo 
mass (solid black line with dotted lines indicating the l-cr 
cluster-to-cluster scatter) to those of L12 and B13. As 
before, at large halo masses, we predict less stellar mass 
for a given halo mass. The discrepancy is particularly 
noticeable around the peak of star formation efficiency, 
Mh 10^^, where the field SHMR is almost 3 times 
what we find in Virgo. Our results here are based on 
the small number (6-7 in total) of very massive galaxies 
in the cluster core, so shot noise introduces considerable 
uncertainty in the abundance matching. Nonetheless, 
our results suggest a lower star formation efficiency for 
the most massive galaxies in the core of Virgo. In the 
conventional picture of hierarchical structure formation, 
the progenitors of these galaxies have always occupied 
high-density regions, so it is plausible that their star for¬ 
mation may have been suppressed even before infall into 
the cluster. 

At lower halo masses, our results extend the trend seen 
between ~ 2 x 10 ^^ and down to halo masses of 

^ 5 X 1 O®M 0 . This appears to conflict with the results of 
B13 at z = 0.1 and z = 1 (thin black and blue curves), 
but only a fraction of the subhalos in the core of the 
cluster have infall redshifts this low, as is apparent from 
Figures IH [TTl or na Instead, our mean results are most 
sensitive to subhalos with z > 1 - 2 . 

Since our abundance matching is based only on Minf, 
independent of Zinf, we do not expect any significant red¬ 
shift dependence in the final results. This is confirmed 
in the top right panel, where we show results for individ¬ 
ual subhalos binned by infall redshift (points with error 
bars; colors correspond to z < 0.5 (blue), 0.5 < z < 1.5 
(cyan), 1.5 < z < 2.5 (green), 2.5 < z < 3.5 (yellow), 
and z > 3.5 (red)). Our derived SHMR is independent 
of Zinf to within the scatter, as expected. The models 
of B13 and L12 predict some evolution in the SHMR, 


particularly at low masses, so this is something we hope 
to test for with the results of the full NGVS stellar mass 
function, as discussed further in appendix [Al 

Uncertainties in the stellar mass function also affect 
the derived SHMR, as illustrated in the lower left panel. 
The four thick (solid) curves show a range of Schechter- 
function slopes and normalizations consistent with the 
pilot region stellar mass function, once photometric er¬ 
rors and population synthesis uncertainties are taken into 
account. Models with shallower slopes and/or higher 
normalizations are in excellent agreement with the re¬ 
sults of B13 at z =0.5-2.5, for halo masses up to 
^ 5 X IO^^Mq, although they still predict a lower SHMR 
at the largest halo masses. The dotted curve shows 
the effect of systematically higher mass-to-light ratios, 
if we have overestimated the importance of bursts in our 
SFHs, as discussed in section |2l This would eliminate 
the discrepancy with the field results, except for the few 
most massive systems (halo masses of l- 2 xlO^^M 0 or 
more), where the cluster-to-cluster scatter is large. 

One shortcoming of our method is the systematic un¬ 
certainty associated with subhalo counting. The bot¬ 
tom right panel shows how the three models presented 
in section lT4l predict dramatically different SHMRs. The 
raw AHF results (model 0) find relatively few surviving 
subhalos in the cluster core, so galaxies are matched to 
objects further down the mass function, producing an 
extremely high SHMR (upper dotted black line). This 
seems inconsistent with the held results of B13 at all red- 
shifts, conhrming that the raw AHF catalog probably un¬ 
derestimates subhalo abundance. The most conservative 
subhalo counting model, model 1 (lower dashed curve), 
predicts a SHMR 5-10 times lower than that of B13 at 
all stellar masses, which also seems implausible given the 
held results. Our intermediate model, model 2 (middle, 
solid curve) matches the results of B13 fairly well where 
the two overlap in mass and redshift, except for the 6-7 
most massive objects, as discussed previously. Assuming 
the normalization in B13 is correct, model 2 therefore 
seems the most plausible method for counting subhalo 
ancestors. We will have the opportunity to test this 
model further when results are available for the entire 
NGVS survey region; at that point the distinct spatial 
and velocity distributions predicted by the three mod¬ 
els (cf. Figure [5]) should provide an alternative way of 
selecting between them. 

Finally, we also consider the effect of an overall mass 
rescaling on the SHMR. The mass of Virgo has been ar¬ 
gued to lie anywhere in the range 0.33-1.2 times our fidu¬ 
cial value 5.76 x 1 O^^M 0 , as discussed in sectional Since 
our simulations were rescaled to the fiducial mass, adopt¬ 
ing a smaller mass would reduce all subhalo masses by 
the same factor. On the other hand, adopting a smaller 
virial radius would increase the fraction of the cluster 
covered by the pilot survey, and thus the fraction of sub¬ 
halos located in this region. Adjusting the mass scalings 
in our simulations, we find these two effects partially can¬ 
cel; adopting a mass 1.2 times larger has little effect on 
the pilot region SIMF, while adopting a mass 0.33 times 
smaller shifts it down by a factor of 2 in subhalo mass. 
The long-dashed curve in the lower right panel of Figure 
|5| shows the effect of this shift; clearly this would also 
help reconcile the results in Virgo with the field SHMR 
determination of B13 at redshifts 0.5-1.5 (thin blue/cyan 
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Figure 8. The derived relationship between stellar mass and halo mass, compared to the results of B13 and L12. The thin colored lines 
show the models of Behroozi et al. for 2 : = 0.1, 1, 2, 3, 4 and 5 (black, blue, cyan, green, yellow, and orange respectively), while the dashed 
lines show the results of Leauthaud et al. for 2 = 0.37 and 0.88 (black and blue lines respectively). The thick black line shows our derived 
SHMR, and the dotted lines show the 1-cr cluster-to-cluster scatter. 


curves), or the measurement of L12 at z = 0.88 (dashed 
blue curve). 

6. CONCLUSION 

Given the complex network of physical processes in¬ 
volved in galaxy formation, simplified descriptions of the 
final outcome, such as abundance matching, provide an 
important middle ground between resolved simulations 
and observational data. In this paper, we have used the 
stellar mass function for the core of the Virgo cluster re¬ 
cently measured by the NGVS to extend the technique 
to unprecedentedly low masses. This exercise is fraught 
with technical challenges. To relate observed galaxies to 
predicted dark matter substructure, we have to deter¬ 
mine how much substructure should survive in the core 
of the cluster (or more precisely, along a line of sight 
through the core), trace this substructure back to the 
point where it last existed as an independent ‘field’ halo, 
and measure the ‘infall’ mass of that halo at that point. 
Abundance matching then relates final stellar masses to 
infall masses, quantifying the net efficiency of star forma¬ 
tion in the halos that were to become cluster substruc¬ 
ture. 

There are uncertainties at each step of this process. 
Substructure can be disrupted artificially in the sim¬ 
ulations by relaxation effects. This problem is rela¬ 
tively easy to correct in an approximate way, by search¬ 


ing through the full merger tree of the cluster to iden¬ 
tify structures that should have survived to the present 
day. On the other hand, this correction introduces sev¬ 
eral uncertainties into the resulting SIMF. The dense 
baryon-dominated components of galaxies will help pre¬ 
serve their surrounding dark matter substructure against 
tidal disruption, for instance, but could also contribute 
to its disruption through strong stellar winds or similar 
feedback effects. Thus the true survival rate for substruc¬ 
ture is not really known. 

A second problem relates to sub-substructure, that 
is to say subhalos of subhalos. Physically, these ob¬ 
jects could correspond to small satellite galaxies of larger 
galaxies, or to individual members of galaxy groups that 
fall into the cluster. In some instances sub-subhalos may 
merge with their parent subhalos before or while the par¬ 
ents fall into the cluster; in other cases they may be dis¬ 
sociated from their parents and survive as independent 
objects. Here again, the precise effect on the final SIMF 
is unclear. 

Finally, we are faced with the usual observational un¬ 
certainties - uncertainties in the global mass and struc¬ 
ture of the Virgo cluster, in the completeness of the 
galaxy detections, particularly at low stellar mass and/or 
low surface brightness, and in the conversion from light 
to mass. These problems have already been addressed in 
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Figure 9. The derived SHMR, compared to the results of B13 and L12. In each panel, the thin solid and dashed lines show the results 
of B13 and L12, with colors and line-styles as in Figure |8] The top left panel shows the mean SHMR derived for model 2, with l-cr 
simulation-to-simulation scatter (thick black line and dotted black lines respectively). The top right panel shows individual mass ratios 
binned by infall redshift (points with error bars; colored as indicated). The bottom left panel shows the effect of alternate fits to the stellar 
mass function (solid lines), as well as the effect of fewer bursts in the SFHs (dotted line). The bottom right panel shows the SHMR derived 
for the two other subhalo counting models, model 0 (the raw AHF results - upper dotted curve) and model 1 (lower dashed curve), as well 
as the effect of reducing the mass of Virgo to 1/3 of the fiducial value (dot-dashed curve). 


Paper I, and will also be discussed further in forthcoming 
NGVS papers. 

Despite these many caveats, if we combine our most re¬ 
alistic model of the SIME with our best estimate of the 
stellar mass function in the core of Virgo, the result is 
a strikingly simple stellar-to-halo-mass relationship, con¬ 
sistent with previous results such as B13 or LI2, but ex¬ 
tending these downwards by almost two decades in halo 
mass or four decades in stellar mass, over the redshift 
range 2 ~ 1-4. Combining our results with those of B13 
and L12, we find evidence for a simple broken power- 
law SHMR over the entire range of stellar mass observed 
outside the Local Group, from IQ^Mq or less to almost 
This two-part fit to the SHMR in turn suggests 
that two main mechanisms may set the net efficiency of 
galaxy formation, one above M* ^ 10"Mq, and one be¬ 
low this mass. 

As with most abundance matching studies, our method 
starts from the assumption that all halos down to some 
mass limit are occupied by galaxies. An alternative pos¬ 
sibility is that galaxy formation is ‘stochastic’, and halos 
of a given mass can contain a wide range of final galaxy 
masses. Methods that combine multiple constraints, 
such as those of B13 or L12, can exclude a large amount 
of ‘stochasticity’ for luminous galaxies. Eor dwarf galax¬ 


ies the situation is less clear, but detailed modeling of the 
Local Group suggests stochasticity may be important 
at lea st in lower-density environments fe.g.. lSawala et ^ 
12014 and references therein). If galaxy formation is very 
stochastic then some halos must be unoccupied, and thus 
the SHMR for the occupied subhalos must be larger, and 
must extend to even smaller halo masses, e.g., down to 
Mh ^ IO^Mq rather than Mh ~ Direct dy¬ 

namical studies of the smallest Virgo dwarfs may help 
constrain this possibility, by placing lower mass limits 
on the dark matter subhalos hosting these galaxies. 

At the lowest stellar masses, we still do not have 
any examples of galaxy formation outside the Local 
Group (nor indeed outside the inner hal o of the Milky 
Way) . The ‘ultra-faint’ dwarf galaxies (IWillm an et al.l 
12004 iBelokurov et al.l 120061 : iZucker et al l 120061) . with 
stellar masses less than IO^Mq, are currently identi¬ 
fied only by star counting, and thus remain invisible 
to us in more distant systems. It is interesting, how¬ 
ever, to consider the implications of our derived SHMR 
for slightly brighter dwarfs in a system like the Lo¬ 
cal Group. At a stellar mass of 3 x 1 O®M 0 , from 
our SHMR we predict an average infall halo mass of 
3 X 1 O®M 0 . Typically, such a system would be stripped 
down to two-thirds of this initial dark matter mass, or 
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2 X 1 O®M0, by the present day. In the Local Group, 
we see fewer than 20 satellites per prim ary above this 
stellar mass limit l)McConnachiel I2012D . whereas from 
simulations, we expect twice as many subhalos over 
2 X 1 O®M0. Thus, we still have a ‘missing satellite’ prob¬ 
lem in the Local Group relative to Virgo: halos that 
would have been occupied by galaxies, had they been 
the ancestors of Virgo cluster subhalos, did not host vis¬ 
ible satellites that survived infall into the Local Group. 
There are several possible solutions to this puzzle. The 
most plausible is that gas cooling and star formation 
are suppressed i n small field halos aft e r reionization 
(lEfstathioul Il992t iBarkana &: LoebI Il999t iBuhock et al.l 

l2000HGnedin fc Kravtsovll2006HSawala et al.ll2014ll : the 

difference in subhalo occupation between the Local 
Group and the Virgo cluster would then relate to the dif¬ 
ferent age distributions of the subhalos’ progenitors. We 
will investigate this possibility further in future work. 

While the results presented here are subject to a num¬ 
ber of uncertainties, the prospects for reducing these un¬ 
certainties in the near future are excellent. First, our 
stellar mass function covers only the core of Virgo, less 
than l/25th of the total area of the NGVS. Thus we 
should expect the uncertainties in the stellar mass func¬ 
tion to decrease dramatically once the full data set is 
analysed. The NGVS also has data at other wavelengths, 
including the near-infrared. We may be able to use this 
to red uce the uncertainties in the stellar mass-to-light 
ratio (jGourteau et al.l [2014( 1. The spatial distribution 
of galaxies across the whole cluster should help identify 
which of our models for the SIMF is correct; any varia¬ 
tion in the spatial distribution with stellar mass will also 
allow us to search for redshift variations in the SHMR, 
as discussed in the Appendix. Currently, redshift infor¬ 
mation for the cluster is extremely limited relative to 
the number of galaxies known; in the cluster core, for 
instance, less than 15% of the galaxies detected in imag¬ 
ing have redshifts in the NASA Fxtragalactic Database 
(NFD). Ongoing work by the NGVS collaboration and 
others (see Paper I, section 6) should help to address 
this, at least for the brighter galaxies. More complete 
spectroscopy would help confirm which of our estimates 
of the SIMF is most realistic, and might allow us to de¬ 
tect redshift dependence in the SHMR. In future work we 
will also complement this additional observational data 


with higher resolution simulations, so that uncertainties 
due to resolution effects are reduced. A detailed study of 
the full galaxy population in Virgo should provide con¬ 
straints on the efficiency of galaxy formation that are 
unique and complementary to any others currently avail¬ 
able. 
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APPENDIX 

ESTIMATING INFALL REDSHIFT FROM SUBHALO ORBITAL PROPERTIES 

The main focus of this paper has been abundance matching in the core of the Virgo cluster, using only the total 
stellar mass function for this region, and ignoring any possible redshift dependence in the SHMR. We expect subhalo 
properties such as infall redshift to vary systematically with cluster-centric position and line of sight velocity, however. 
Thus results from the full NGVS survey and/or additional spectroscopic data should allow us to construct multiple 
samples with different mean Zi^f, and thus refine our understanding of the correspondence between galaxies and 
substructure. We explore this possibility further in this appendix. 

The correlations between subhalo locati on, orbital properties and infall redshift have been noted previously in many 
semi-analytic ie.g.. lTavlor fc Babnlll2005^ and numerical (e.g., Gao et al. 2011, 2012; Contini et al. 2012 and references 
therein) studies. Figure ITOl shows the correlation between radial position at z = 0 (left panel), or infall redshift (right 
panel), and normalized orbital energy at ^ = 0. Subhalos that merge at early times generally occupy lower energy 
orbits, and the lowest energy orbits are found at smaller radii. The converse is not quite true, however - there is a 
population of objects on weakly bound orbits that are spread over a range of radii (top part of the distribution in the 
left panel). Further examination shows that these are systems on their first passage in to pericenter, so they are recent 
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Figure 10. Subhalo orbital energy Eo^hi normalized to the square of the main halo’s circular velocity at the virial radius, Uc, versus (3D) 
radius r normalized to the virial radius rvir (left) and infall redshift 2;inf (right). For clarity, only one subhalo in 20 is plotted. Red data 
points and error bars indicate the mean and 1-a scatter in evenly spaced bins. 



Figure 11. Intrinsic properties in terms of observables. Left: normalized orbital energy as a function of projected separation 

Rproj and velocity offset (normalized to the virial radius and mass of the main halo respectively). Right: infall redshift in the 

same phase space. For clarity, only one subhalo in 100 is plotted. The vertical dashed line in each panel shows the approximate extent of 
the pilot region. 


mergers. 

Orbital energy is not directly observable, however; instead we must look for patterns in projected radius -Rproj and 
line of sight velocity v\os- Figure ITT] shows the distribution of subhalos in the space of observable quantities, colored 
by relative orbital energy (left panel) or infall redshift (right panel). Here we see the separation into recently merged 
and older material familiar from studies of the phase-s pace distribution of cluster galaxies fe.g.. iMahaian et al.ll201l1 : 
iHaines et al.ll20l3 lOman et al.ll2013l : iHou et aT1l2014l l. Orbital energy is particularly well correlated with position in 
phase space. The distribution of infall redshifts is more mixed, but we may still be able to extract useful information 
from the observables. 

The trend of decreasing Zint with radius is well-known, and has been used for decades to study galaxy transformation 
in clusters. Data from the central region of NGVS, or similar studies of cluster cores, cannot probe these large-scale 
radial gradients. There should also be a trend with velocity offset from the cluster center, however, particularly at 
small projected radii. If we define samples of subhalos by their line of sight velocity offset from the mean cluster 
velocity uios, we can obtain different distributions of mean orbital energy and Figure [T^] illustrates this for cuts at 

uios/t'esc = 0.2 and 0.4, where Vesc = ^/2\(j){x)\ is the local escape speed of the cluster. We see that the high-velocity 
sample in particular contains a large excess of loosely bound systems with z-mi < 0.5. 

Thus with coverage of the entire projected area of the cluster, and/or more complete spectroscopic information, we 
will be able to separate Virgo galaxies into subsamples with different mean infall redshifts. This should increase the 
sensitivity of abundance matching to redshift variations, reduce the uncertainties in the modeling of the SIMF, and 
constrain stochasticity in subhalo occupation. 
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